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In this paper we investigate the properties of self induced transparency (SIT) solitons, propagating 
in a A-type medium. It was found that the interaction between SIT solitons can lead to trapping 
with their phase preserved in the ground state coherence of the medium. These phases can be 
altered in a systematic way by the application of appropriate light fields, such as additional SIT 
solitons. Furthermore, multiple independent SIT solitons can be made to propagate as bi-solitons 
through their mutual interaction with a separate light field. Finally, we demonstrate that control 
of the SIT soliton phase can be used to implement an optical exclusive-or gate. 

PACS numbers: 42.65.Tg, 42.79.Ta, 42.50.Gy 



I. INTRODUCTION 

Strong and coherent interactions between A-type 
three-level atomic media and resonant light pulses have 
been a topic of much interest recently. Driving this in- 
terest are various optical communications applications, 
such as optical buffers and phase shifters for optical 
routers [U El • Of particular interest is all-optical switch- 
ing based on solitons to realize optical gates ([3 and refer- 
ences therein) . Generally, the Kerr effect is chosen as the 
nonlinear interaction between solitons due to its presence 
in many different media. However, the low value of typ- 
ical Kerr coefficients implies that intense optical pulses 
and/or long interaction lengths must be used. In con- 
trast, self-induced transparence [4 (SIT) provides soli- 
tons that are based on a strong, resonant optical non- 
linear effect, which could, in principle, reduce interac- 
tion lengths and pulse intensities required for continuous 
variable optical computing. This property makes it an 
attractive nonlinear optical medium for optical comput- 
ing. 

The interaction of a single SIT soliton with a con- 
trol pulse has been investigated by Vemuri et al. They 
showed that a single SIT soliton can be cloned in a A-type 
medium by applying a weak control light field resonant 
with the 1 2)- 1 3) transition (see Fig. [T]) [5j. In doing so, 
the spatial region of the medium where the SIT soliton 
and control pulse overlap is left with a significant popu- 
lation in 1 2) and a ground state coherence extending to 
either side of the region [6, 7 , which we refer to as the 
trapped soliton. 

The trapped soliton presents an intriguing possibility: 
namely, that linear and nonlinear optical operations may 
be performed on the trapped soliton using resonant light 
fields. Normally, operations between different solitons 
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must be carefully timed so that the overlap in the lin- 
ear (e.g., beam-splitters) or nonlinear optical elements is 
maximized. However, in the case of a trapped soliton, the 
timing problem is much less critical and computational 
operations to be performed on the stationary coherence, 
rather than on the traveling soliton. 

Unfortunately, very little is known about the quantum 
or classical interactions between trapped solitons, prop- 
agating SIT solitons, or other light fields. Recent work 
has shown that the presence of a trapped soliton is suf- 
ficient to capture additional solitons [HIH], however, the 
mechanism of such trapping, and the interaction of the 
trapped solitons with light fields resonant with the |2)-|3) 
transition remains unexplained. In what follows, we go 
beyond these results by examining the influence of addi- 
tional control pulses and SIT solitons on the state of the 
medium and the phase of the emitted light. 

The details of the capture process are explained and 
we highlight this by examining three representative cases. 
The interaction of a control light field with both single 
trapped solitons and multiple solitons is examined. It 
was found that multiple trapped solitons can behave like 
bi-solitons. Further, we show that the capture process 
and interactions with control light fields manipulates the 
phases of the trapped soliton. Finally, by choosing the 
soliton phase as an information carrier, we demonstrate 
a series of operations that implement an all-optical XOR 
gate. 



II. THEORY 

To study the propagation and interaction of light in a 
three-level medium, where each atom has energy states 
as shown in Fig. [l] we employ a density matrix ap- 
proach in the rotating wave approximation (RWA) [9]. 
The medium's constituent atoms have one excited state 
and two ground states. The |1)-|3) and |2)-|3) are dipole 
allowed radiative transitions, while the |1)-|2) transition 
is dipole forbidden. The medium is discretized along the 
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FIG. 1: A-three level system with two ground states and one 
excited state. A transition from state |1) to |3) can be per- 
formed with a SIT soliton. The transition from |2) to |3) can 
be performed with control light field. A radiative transition 
between |1) and |2) is dipole forbidden. 



propagation axis, creating a 1-D grid. At each grid lo- 
cation along the propagation axis, the medium's state 
is represented by an independent density matrix. This 
leads to the following set of differential equations for each 
grid point, the solution of which provides the time and 
space-dependent population densities and coherences be- 
tween populations for the whole medium. 

Pu = - 1 (Xi3Pi3 - xtaPis) + PssRsi (1) 

P22 = - 1 (X23P23 - XhPh) + P3SRS2 (2) 

P33 = t (X13P13 + X23P23 - X^3Pl3 " X53P23) 

-(i^31+i^32)p33 (3) 

Pl2 = \ (xt3P53 - X23P13) (4) 

Pl3 = \ (Xl3 (P33 - Pll) - X53P12) - /^13Pl3 (5) 

P23 = I (X23 (P33 - P22) - X\zP\^ - P22.P22. (6) 

The diagonal elements of the density matrix, pa^ are 
the population densities, normalized such that Tr(p) = 1, 
of the energetic states of the medium. The off-diagonal 
elements, pij^ are the coherence terms that are propor- 
tional to the polarization of the medium. The decay 
constants R^^i and i?32 account for spontaneous emis- 
sion from 1 3) to |1) and |2), respectively. Decay of the 
coherence between states i and j is given by the rate, 
Pij^ representing non-radiative processes, such as colli- 
sional dephasing. The relationship between the Rabi fre- 
quencies, Xiji arid the applied light pulses with a field 
envelope, is given by 
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where pis and /i23 are the transition dipole moments of 
the two radiative transitions. 

The light fields, which are assumed to be resonant with 
their respective transitions. The duration of the light 
pulses are assumed to be much shorter than the spon- 
taneous emission rates, but long enough for the slowly 
varying envelope approximation (SVEA) to hold. We 



also require that the life-time of |3) is sufficient to al- 
low stimulated emission processes dominate spontaneous 
emission. The temporal profile of the light field envelopes 
are given by sech-squared functions, as SIT solitons as- 
sume this form naturally [4 . 

Applying the SVEA and limiting ourselves to one spa- 
tial dimension (plane wave propagation into the positive 
z-direction), the Maxwell equations are given by: 



—- = Z — Np23p23 
CtZ Cq 



(9) 
(10) 



with ks and kc the wavenumbers of the SIT and coupling 
fields, N the density of the active medium and eo the 
electric permittivity of the vacuum. 

The coupled set of density and Maxwell equations 
Eqs. [l][6j |9] and [lOj were solved numerically using both 
a fourth order Runge-Kutta algorithm. All calculations 
begin with the |1) state full occupied (pn = 1) and no 
coherence between any of the levels. 

in our calculations we use a value for the product of 
Np = 1.7 • 10~^^ Cm~^ [10 , which corresponds to prob- 
ing the transition from state 4^P to state 3^D of calcium 
in a gas-cell with a 30 mBar partial pressure of calcium. 

After a number of solitons and control pulses have in- 
teracted with the medium, we examine the spatial struc- 
ture of the ground state coherence (P12), the phase of the 
ground state coherence, and the phase of light emitted by 
the medium. In general the phase of the light envelope is 
given by the argument of Ei = Aiexp{i(j)i) . Likewise the 
phase of the ground state coherence is given by the ar- 
gument of pi2 = Apexp{i6). The phases are defined with 
respect to an initial control pulse, which we consider to 
have a phase of zero. In other words, the control pulse 
has a purely real field envelope. 



III. RESULTS 

We begin by partially repeating the work of Vemuri 
et al. [5 , where a single SIT soliton is trapped by the 
application of a single control pulse. This occurs because 
as the SIT soliton excites population into |3), the control 
pulse stimulates emission into |2). The result is that 
the control pulse is amplified at the expense of the SIT 
soliton, which is absorbed by the medium. As shown by 
Loiko et al. [7 , in this process, the medium has stored 
the sech-squared shape of the SIT soliton in the medium 
as a sech-squared shaped distribution of population in 
the second ground state. We observe that, in addition to 
the shape preservation, the phase of the stored soliton is 
also stored in the ground state coherence of the medium 
(see red trace in Fig. |2|. More precisely: 



^13 



(11) 



where ^13 is the phase of the SIT soliton. This re- 
lationship comes from the fact that the coherence of 



3 



q: 0.0 



-0.5 




100 200 300 400 500 

Position (jjm) 



FIG. 2: State of the ground state coherence after the storage 
of a single SIT sohton by a small coupling pulse (red line). 
Ground state coherence after the application of a 7r-area con- 
trol pulse to a medium containing a stored SIT soliton (blue 
dotted line). The excited region is shifted deeper into the 
medium and the ground state coherence undergoes a tt phase 
shift. 



the medium is proportional to the polarization of the 
medium, which opposes the light fields. Our calcula- 
tions show another phenomenon: The blue dotted trace 
in Fig. [2] shows the resulting ground state coherence af- 
ter a TT-area control pulse has interacted with the trapped 
SIT soliton. Note that the 7r-area control pulse shifts the 
SIT soliton's phase by tt, while also shifting the location 
of the trapped SIT soliton. However, control pulses with 
an area less than tt shift the location of the ground state 
coherence without altering its phase [llj. 



A. Interactions between SIT solitons 

We now turn to the case where a second SIT soliton 
is injected into the medium after the capture of a first 
SIT soliton. In this case, the second SIT soliton encoun- 
ters a medium with populations and coherences described 
above and illustrated in Fig. |2j As will be described in 
more detail, we observed that, in general, these addi- 
tional SIT solitons are captured as well upon "collision" 
with the excited region, and that this leads to the cre- 
ation of additional excited regions. An important feature 
of this process is that the second and subsequent SIT soli- 
tons are captured without a control pulse being injected. 

This can be understood by examining the behavior of 
the SIT soliton and the A-system, illustrated in Fig. |3] 
Fig. [sja) shows the intensity of the SIT solitons as a 
function of time and space, while Fig. [sj^b) shows the 
intensity of control light fields. An initial SIT soliton is 
injected into the medium at time A, which travels slowly 
through the medium. A short time later, a control pulse 
is injected, it is however, too weak to be visible on the 



scale of Fig. |3jb). The control pulse overtakes the SIT 
soliton at time/position where it is amplified at the 
expense of the SIT soliton, producing the excited region 
illustrated in Fig. [2] Note that the amplified control pulse 
exits the medium without further perturbation. 

A second SIT soliton is injected into the medium at 
time C, however, no control light fields are injected. Nev- 
ertheless, the SIT soliton is trapped in the medium by 
the collision with the excited region at time/position D. 
When the SIT pulse enters the first half of the excited 
region, the ground state coherence ensures that there is 
simultaneous stimulated emission on both the |1)-|3) and 
1 2)- 1 3) transitions. However, the light on the control tran- 
sition is out of phase with the original control pulse, in- 
dicated by a negative value in Fig. [sj^b) at location D. 
This control pulse is amplified at the expense of the SIT 
soliton, creating a new excited region in the medium. 

At location the control light field enters the second 
half of the excited region, where the sign of the coherence 
is the reverse to the first half. The population transfer 
due to the control pulse destroys this coherence, generat- 
ing a new SIT soliton, which we refer to as a transitory 
SIT soliton. In doing so, the control pulse is completely 
absorbed and a new control light field with the phase of 
the original control pulse is emitted by the residual pop- 
ulation and ground state coherence. This control pulse 
captures the newly generated SIT pulse, being ampli- 
fied in the process and exits the medium (location F). 
The second capture process generates another excitation 
region. As a result, the medium now has a differently 
shaped excitation region, which is wider than before. 

It is also important to note that the soliton collision 
results in a control pulse exiting the medium. This is 
a general result, and also occurs if subsequent collisions 
occur. 

As described above, a transitory SIT soliton is emitted 
during the collision and subsequently recaptured by con- 
trol light fields that are also emitted during the collision 
process. Thus, one would expect that the coherence gen- 
erated during the capture of the transitory SIT soliton 
to replicate that of the second SIT soliton capture. How- 
ever, the original coherence, which is the source for all the 
control light fields after the first SIT soliton is captured, 
is an odd function with respect to the center location of 
the captured SIT soliton. Thus, the control light fields, 
emitted from each side of this location during the colli- 
sion, are out of phase with each other. Consequently, 
these two light fields destructively interfere with each 
other. In this case, the two SIT solitons are identical, 
leading to a near complete destructive interference of the 
control light field. The control light that subsequently 
captures the transitory SIT soliton is generated from the 
residual excited population and coherence, which has a 
phase identical to the original control pulse, thus, the 
coherence associated with the transitory SIT soliton is 
identical to that of the first SIT soliton. 

This is highlighted by examining three cases: Ppr = 
Ptr, Ppr > Ptr^ Ppr < Ptr^ whcrc Ptr IS the peak intensity 
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FIG. 3: Electric field amplitude of light resonant with the |1)- 
|3) transition, as either an injected SIT soliton, or emission 
from the medium, (a) Electric field amplitude (V.m~^) of 
light resonant with the |1)-|3) transition and (b) Electric field 
amplitude (V.m""*^) of light resonant with the |2)-|3) transi- 
tion, as either an injected control pulse, or emission from the 
medium. Event A: the first SIT soliton is introduced to the 
medium. Event B: the SIT soliton and control pulse over- 
lap, resulting in the SIT soliton's capture and the coupling 
pulse's amplification. Event C: the second SIT soliton is in- 
troduced to the medium. Event D: the SIT soliton collides 
with the first half of the excited region, emitting a negative 
control pulse. Event E: the negative control pulse captures 
the SIT soliton and is amplified. It then destructively inter- 
feres with the ground state coherence of the second half of the 
excited region, resulting in the emission of a transitory SIT 
soliton. Event F: the remaining control light field is amplified 
by the transitory SIT soliton, capturing it in the process and 
extending the original excitation region. 



of the SIT soliton that is trapped by a control light field 
and Ppr is the intensity of the SIT soliton that collides 
with the ground state coherence, which was written into 
the medium due to the capture of Ptr- 

CclSG ll Ppr ~ Pfr 

For the case of equal intensity SIT solitons, the sec- 
ond SIT soliton collides with the excited region and, as 
a result, does not penetrate as far into the medium as 
the first SIT soliton (see Fig. [4| in comparison to Fig. |2j 
The act of colliding emits a transitory SIT soliton, which 
expands the excited region deeper into the medium by 
approximately twice its original full width half maximum 
(FWHM) width. The two peaks in p22 are separated by 
4 times the FWHM width of the original excited state re- 
gions. Note that P22 falls to zero between the two peaks, 
however, as we will show later, these two regions remain 
coupled. 

Note also that the phase of the coherence appears to be 
inverted with respect to the original capture. The inver- 
sion occurs because the control light field that captures 
the second SIT soliton is generated by the coherence on 
the left hand side of the excited region in Fig.|2|. This 
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FIG. 4: p22 (red trace) and ground state coherence (blue 
trace) through out the medium after two SIT solitons with 
equal peak intensity are captured. 



light is out of phase with respect to the original control 
pulse (we also refer to this as a negative control pulse), 
thus, from eqn. [nj a shift of TT is expected. 
Case 2: Ppr < Ptr 

As in the previous case and shown in Fig. [5j the second 
SIT soliton is indeed captured, however, the resulting 
state of the medium is somewhat different. The lower 
amplitude and slower traveling second SIT soliton emits 
a lower amplitude and slower traveling transitory SIT 
soliton during the collision. Consequently, the second 
P22 peak is not shifted as deeply into the medium (60% 
of its FWHM width) compared to case 1. The distance 
between both peaks is about 6.5 times the FWHM of the 
original excited region, which is greater than for case 1. 
The greater separation is due to the fact that the second 
injected SIT soliton travels slower, and is thus, captured 
over a shorter distance by the emitted control light fields. 

As before, the excited population is surrounded by 
ground state coherence with the ground state coherence 
undergoing a tt phase transition when P22 returns to zero 
between the two peaks. This indicates that, as in the 
previous case, the transitory SIT soliton is captured by 
a positive coupling pulse. Given case 1, where equal am- 
plitude SIT solitons result in a positive coupling pulse, 
this is expected. The relatively less intense control fields, 
emitted during the capture of the second SIT soliton on 
the left hand side of the original capture location, de- 
structively interfere with the more intense control light 
fields emitted from the right hand side of the capture lo- 
cation. The end result is a control light field with the 
same phase as the original control pulse. 

CclSe 3 1 Ppr ^ Pfr 

In this case, the second SIT soliton is more intense 
and has a shorter temporal duration than the prepara- 
tory pulse, resulting in quite different behavior compared 
to the previous two cases. The second SIT soliton is 
still captured, however, it generates a fast moving, high 
amplitude, transitory SIT soliton that is captured much 
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FIG. 5: p22 (red trace) and ground state coherence (blue 
trace) through out the medium after two SIT sohtons are 
captured. In this case, the peak intensity of the second SIT 
soUton is less than that of the first. 



FIG. 6: p22 (red trace) and ground state coherence (blue 
trace) through out the medium after two SIT solitons are 
captured. In this case, the peak intensity of the first SIT 
soliton is less than that of the second. 



deeper in the medium than the location of the original 
excited region (see Fig. 6|. Furthermore, the intensity of 
the negative control pulse, generated during the capture 
of the second SIT soliton from the left hand side of the 
original coherence, is large enough to overwhelm the pos- 
itive control light field, emitted from the right hand side 
of the original excited region. As a result, the transitory 
SIT soliton is captured by a negative control pulse and 
the phase of the ground state coherence remains constant 
between the two P22 peaks. 

Because the transitory SIT soliton has a higher inten- 
sity than the first SIT soliton, it appears as though the 
second SIT soliton travels through the excitation region 
associated with the first SIT soliton without any inter- 
action, but, this is not the case. The second SIT soliton 
is captured shallower in the medium than both the pre- 
ceding cases, due to the strong control field amplification 
associated with the high peak intensity of the SIT soli- 
ton. However, the emitted control light field and transi- 
tory SIT soliton removes much of the excited population 
associated with the capture, leaving the first P22 region 
smaller than expected. As a result of these interactions, 
the second SIT soliton is captured earlier, while the dis- 
tance between the two excited regions is 14 times the 
FWHM of the original excitation region. 



B. Interaction between multiple trapped solitons 
and a control light field 

As shown in Fig. [2j injecting a control light field alters 
the location of the trapped soliton. In this section, we 
investigate this behavior after one SIT soliton has been 
trapped by a control pulse and up to four subsequent SIT 
solitons have been captured by the pre-existing ground 
state coherence. A control light field is then applied and 
the behavior of the trapped solitons is observed. In Fig.[7| 



the electric field, Es^ throughout the medium is shown 
for several different situations. 

Fig. [7|^a) shows two SIT solitons propagating in re- 
sponse to the application Ec. The first SIT soliton is 
launched into the medium at tq = 0.1 ns and is captured 
by a small control pulse. At tq = 0.4 ns a second SIT 
soliton is launched into the medium, colliding with the 
ground state coherence pi2 structure (lower) as left over 
from the first SIT soliton. From t = 0.8 ns on a continu- 
ous wave control light field Ec = 0.5 Es is launched into 
the medium. This light field interacts with the ground 
state coherence structure as left by the two SIT solitons 
in the medium and generates a light field Es that prop- 
agates as a pair of coupled solitons in the positive z di- 
rection. The ground state coherence (not shown in the 
figure) moves together with the pair as well. 

More than two SIT solitons can be stopped and stag- 
gered in the medium, which is shown in Fig. [7| where 
(b) shows the case of three SIT solitons, (c) four and (d) 
five solitons stored. As can be seen, when the cw control 
field is injected from t = 0.8 ns, only pairs of solitons 
are formed and that two of these pairs do not cross to 
form larger structures. Furthermore, for an odd number 
of SIT pulses, one of the pulses (the last as seen from the 
entrance of the medium) does not pair up. 

This behavior can be understood as follows. The con- 
trol field, which is absorbed by the first region of excited 
ground state coherence, is amplified by the second re- 
gion and, as a result, has approximately the same am- 
plitude after passing the excited regions (consisting of 
two trapped pulses) as before entering it. Furthermore, 
where multiple SIT pulses were trapped in the medium, 
each new pulse creates an excited region of ground state 
coherence that is shorter than the previous one. The con- 
trol pulse interacts with these excited regions and creates 
a field Es, where the power in the pulse is proportional 
to the total amount of ground state coherence with which 
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FIG. 7: The electric field Es (V.m~^) throughout the medium 
in a set of experiments where two (a), three (b), four (c), 
and five (d) SIT solitons are trapped in the medium. Sub- 
sequently, from t=0.8 ns, a continuous wave control field is 
launched into the medium that interacts with the ground state 
coherence structure to form sets of bi-solitons and single soli- 
tons which propagate in the +2;-direction. 



the control pulse has interacted. The propagation speed 
of these created pulses is proportional to the power in 
the pulse. Therefore, the second soliton pair (as counted 
from the entrance of the medium) propagates faster than 
the first bi-soliton because it has a larger excited region 
and, as a result, generates more optical power. When an 
odd number of SIT pulses are inserted, the control field 
creates a single pulse in the Es field that propagates 
faster than the closest soliton pair. 

Other systems (with different physical processes at 
work) in which such similar soliton structures propagate 
have been reported by Maruta and Hause [121 US] where 
these soliton structures were named bi-solitons. How- 
ever, such bi-solitons in a A-type system have, to our 
knowledge, not been reported. 

C. Optical logic gates 

Here, we show that an all-optical XOR gate can be im- 
plemented using combinations of SIT solitons and control 
pulses. In the continuous variable regime, using the SIT 
soliton phase as the information carrier, a XOR operation 
is given by |(/)i,(/)2) + (/>2)- The logic operation 

begins with the cloning and capturing of the first SIT 
soliton (01 ), by the application of a control pulse. The 
phase of this SIT soliton is transferred to the ground state 
coherence of the medium at the cloning location, storing 
it. Thereafter, the phase of the coherence is inverted by 



the application of a 7r-area control pulse. Finally, the 
SIT soliton with the phase to be operated on (02 ) is in- 
serted into the medium. This SIT soliton collides with 
the stored coherence, generating a control pulse possess- 
ing the desired sum-phase 0i +02- 

Note that, in addition to encoding the result of the 
XOR operation in the phase of the emitted control pulse, 
the ground state coherence of the medium also stores the 
result of the operation. The stationary nature of the 
ground-state coherence then allows multiple operations 
to be performed. Furthermore, because each operation 
results in the medium emitting a pulse of light that en- 
codes the result of that operation, conditional operations 
become feasible. 

It is interesting to estimate the number of optical logic 
operations that might be performed on a trapped SIT 
soliton. The limit to the number of operations that can 
be performed on a stored soliton can be estimated by 
comparing laser pulse repetition frequencies and ground 
state coherence times. An ultracold atomic ensemble 
quantum memory gate with a 1 ms decoherence time was 
recently demonstrated [14 , while laser repetition rates 
exceeding 300 GHz have also been demonstrated [15]. 
This provides an upper bound of 300 million operations 
on a stored SIT pulse. 



IV. CONCLUSIONS 

We have analyzed the behavior of SIT solitons in a 
three-level medium by solving the density matrix and 
Maxwell equations. We have found that multiple SIT 
solitons can be stored as a ground state coherence in 
the three level medium by either the application of a 
control light field or by colliding the SIT soliton with a 
pre-existing ground-state coherence. We have analyzed 
how the phase of the ground state coherence is related to 
the phase of the SIT soliton and the control pulse. The 
process of capturing SIT solitons has been investigated 
and has been shown to be the result of a complex process 
of repeated SIT soliton absorption and emission. 

The propagation of trapped SIT solitons, driven by a 
control light field, has been investigated. As a result, 
bi-soliton behavior were observed for the first time in 
SIT solitons. Finally, we have shown that the phase of 
the SIT soliton can be used to encode information and 
demonstrate the operations required to produce an all- 
optical XOR gate. We note that the results from every 
operation on the ground-state coherence is stored in the 
ground-state coherence and emitted as a control light 
field pulse. This provides a convenient mechanism for 
implementing conditional operations. 
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